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PREFACE 


In LACIP 'hare III, one of the largest sources of proportion estimation error 
for a ,r interest in Land Satellite (Landsat) data was boundary pixels 

and the associated method for handling them. These are the pixels partly 
within the class of interest and partly not because of the coarseness of the 
resolution of the Landsat sensor. This memorandum describes many of the 
methods and approaches for dealing with the influence of boundary pixels on 
proportion estimation, A general framework is presented for viewing the 
problem, and the results of a small study are presented. The study of the 
effect of boundary pixels in proportion estimation is very new; almost all 
the methods presented are untried. They represent the most current thoughts 
for attacking the problem. 

The principal author and project coordinator, D. T. Register, originated the 
statistical sampling approach that provides the general framework. A. L. Ofia 
conducted most of the analysis of the small study in North Dakota and Kansas. 
Analysts B. B. Schroder, B. A. Tolbert, C. W. Haynes, and C. L. Dailey as- 
sisted in the study by examining thousands of pixels to decide whether they 
were boundary or not. 
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1. INTRODUCTION 


A detailed examination of the segment Procedure 1 proportion estimation tech- 
nique that was used in Phase III of the Large Area Crop Inventory Experiment 
(LACIE)(ref. 1) has shown that the largest single source for small -grain pro- 
portion error is the analyst's misidentification of picture elements (pixels) 
(ref. 2). Sometimes small grains were thought to be nonsmall grains (omission 
error), and at other times nonsmall grains were identified as small grains 
(commission error). The omission error rates were much higher than the com- 
mission rates, this resulted in an underestimation bias. The three major 
causes of analyst labeling errors were found to be (in order of importance) 
boundary pixels, abnormal signatures, and inadequate acquisitions. Approxi- 
mately 40 percent of all mislabeled pixels were associated with boundaries. 
Reference 2 defined a pixel as boundary if it was spatially located so as 
to be only partly within a small-grain field;^ that is, if it lay on the 
perimeter of a small-grain field. The intent of this memorandum is to 

1. Report on the past importance of boundary pixels in Phase III and anti- 
cipate tneir importance for the future. 

2. Discuss some possible approaches for handling boundary pixels to reduce 
the bias toward underestimation of the crop of interest. 

3. Present the results of some experimentation into one of these approaches. 

4. Make recommendations on what considerations should be used in planning 
a boundary pixel research program. 


In the case of strip/fallow areas, the pixels were divided into two kinds - 
those whose strips were wide enough to be discernible in the Landsat imagery 
and those whose strips were so narrow that no stripping was discernible. If 
the strips were not discernible, the pixels were not considered as boundary 
nor in error. 
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2. IMPORTANCE OF BOUNDARY PIXELS 


The singlemost important cause for the mislabeling of pixels in Phase III 
of LACIE was the boundary pixel. Table 1 presents the results of a detailed 
dot-by-dot examination of the labeling of type 2 dots. All the results are for 
Phase III of LACIE. The scope of the study included all the intensive test 
sites (ITSs) and a sample of blind sites drawn from each of five sele'ted 
states. The results show that 66 percent more boundary pixels were Icibeled 
nonsmall grains than were labeled small grains. The percentages of error shown 
in table 1 do not average to the 40 percent mentioned in section 1. The 
percentages in the table refer to errors for which no other cause than bound- 
ary could be given. The 40 percent included all errors associated with 
boundaries. There were in some cases multiple causes. 

TABLE 1.- DISTRIBUTION OF BOUNDARY PIXELS AND 
THEIR CONTRIBUTION TO LABELING ERROR 


Dau set 

No. of 
segnents 

No. of type 2 
boundary 
dots 

No. boundary 
dots 

analysts labeled 
small grains 

No. boundary 
dots 

analysts labeled 
nonsamll grains 

Percentage at 
all pixels 
that are 
boundary 

Percentage of 
oailsslon 
errors due to 
boundaries 

Percentage of 
coaaltslon 
errors due to 
boundaries 

Minttr ITSs 

13 

112 

39 

73 

12.4 

21.4 

8.6 

Spring ITSs 

7 

($ 

30 

35 

12.4 

12.9 

8.8 

HI nnesota 

6 

S« 

20 

36 

16.0 

40.6 

44. < 

Montana* 

10 

55 

28 

37 

6.9 

21.1 

29.4 

North Dakota* 

IS 

120 

43 

77 

11.7 

28. 9 

22.1 

Oklahoma 

12 

80 

27 

53 

10.3 

22.1 

14.3 

Colorado* 

6 

29 

11 

18 

7.3 

25.0 

100.0 


•Th€$* « 9 n»nt$ conUlwd sone itrtp/f«llow fields which were so narrow thjt they could not he discerned In the 
Imeqery. Pixels In these fields were not considered as boundary or In error. 
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The effect of boundary pixels can be expected to grow In Importance In the 
Agricultural and Resources Inventory Surveys Through Aerospace Remote Sensing 
(AgRISTARS) Project as interest is directed to inventories of multiple crops 
in areas with smaller fields. In particular, for corn and soybeans, there 
will be tliree kinds of boundaries: corn/other, soybeans/other, and corn/soy- 

beans, whereas in the past, only small grain/other boundaries were considered 
In addition, the sizes of the fields in the corn and soybean areas are gen- 
erally smaller than those for wheat in the U.S. Great Plains. The following 
table illustrates how the number of boundary pixels could increase as more 
crops are estimated in areas with smaller fields. The table lists five sites 
within the corn and soybean region; for each segment, the percentage of in- 
terior pixels is given. For this table to demonstrate the effect o^ more 
crops, interior pixels are those which are wholly within the same category 
on each acquisition. Each crop was considered a separate category. 


Segment 

State 

Percentage of pixels 
interior to fields 

146 

Kentucky 

45.0 

812 

Missou'i 

46.4 

824 

11 1 inois 

44.0 

883 

Iowa 

42.1 

886 

Idaho 

52.2 


2 

These five segments were selected because of availability of data and are 
only for illustration. They may not be completely representative of the 
corn/soybean area. 
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3. APPROPRIATE DEFINITIONS FOR BOUNDARY PIXELS 

Conmonly, the term boundary pixel refers to a pixel which lies on the spatial 
demarcation between two categories of Interest on a single Landsat acquisi- 
tion. Though this definition Is appropriate for visualization, It Is not 
a working definition. A working definition must specify the method for 
deciding whether or not a pixel Is a boundary. This point Is Important 
because the designation of pixels as boundary can and w H vary considerably 
depending upon the method utilized. Among the possible methods which can be 
chosen are clustering algorithms, classification algorithms, analysts, and 
aircraft photographs with ground truth. Because the aircraft photographs 
are not routinely available for use In estimation of crop segment proportions, 
they will not be further considered In this memorandum as a method of deter- 
mining which pixels of a scene are boundary. 

There are basically two issues In regard to which pixels are boundary. The 
first issue can be summarized by the question: How much of a boundary does 
a pixel need to be before It Is labeled "boundary"? The second, maybe not 
so obvious. Issue is the following: In which Landsat acquisition Is the 

pixel a boundary? Because of the misregistration between Landsat passes 
over a segment, a given pixel can migrate back and forth between and on the 
boundary of two fields. Thus, an adequate answer to whether or not a pixel 
is boundary requires specification of a particular acquisition. Sometimes, 
a set of acquisitions Is specified. In the Label Identification from Statis- 
tical Tabulation (LIST) approach (ref. 3), a set of four acquisitions Is 
specified. If a pixel migrates back and forth or is on the boundary of two 
fields as determined by an analyst in these acquisitions, then the pixel Is 
designated a boundary. In Procedure 1, a similar approach Is taken: The 
analyst is asked to label for type 1 dots only those pixels which remain 
within the same field over the four (or fewer) acquisitions chosen for proc- 
essing. The implication for type 1 dots Is that all pixels not labeled are 
boundary. 

^In procedure 1, the type 1 dots are first used for starting a clustering 
algorithm and then for labeling the resulting clusters by the nearest- 
neighbor rule. With type 2 dots, all are identified and used for 
stratified areal estimation. 
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Two proposals have been made for use by analysts in this determination. The 
first proposal is that a pixel should be called a boundary pixel if an analyst 
has difficulty ascertaining to which of the two contiguous categories on which 
the pixel lies he should assign the pixel. The second proposal is that a pixel 
should be called boundary if the analyst can determine that any part of the 
pixel is on the spatial boundary of two categories. The second proposal is 
more liberal than the first and will result in a larger number of boundary 
pixels. One must keep in mind here, however, that the importance of the num- 
ber of boundary pixels is secondary to that of consistency in determination. 

Several analysts were asked which of the two proposed working definitions 
would provide the most consistent responses for denoting pixels of a segment 
as boundary. The overwhelming opinion was that the liberal definition would 
provide more consistent responses. 

The liberal definition is the one adopted for use in the study described in 
section 4.3.2. As for machine algorithms, such as clustering and classi- 
fication, they would objectively and consistently apply the rules for deter- 
mining boundaries that they are programmed to follow. Other procedures may 
require their own particular definitions, such as the LIST approach. In the 
LIST approach, whether or not a pixel is boundary is determined by an analyst's 
response to a given set of questions. In summary, the liberal definition 
appears to be an appropriate definition for an analyst to use. However, 
special procedures or algorithms will require their own definitions. All of 
these can be expected to give differing results. In addition, the designa- 
tion boundary for a pixel is only applicable to a particular acquisition or 
set of acquisitions. 


5 


4. TECHNIQUES FOR REMOVING THE ESTIMATION 
BIAS DUE TO BOUNDARY PIXELS 
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Four bask types of techniques have been suqqested for handling boundary 
pixels: cluster-based technique, maxiinum-l ikel ihood-based technique, 

statistical s.implinq technique, and multi category labeling. 

^ 1 CLUSTER-n ASED TEC HNIQU ES 

The cluster-based terhniques diverge into two separate approaches. The dis- 
tinctive difference between the two is the means of determining which pixels 
are boundary. The means may be either automatic (e.g., allowing a clustering 
algorithm to decide which pixels are boundary) or manual (e.g., in the case 
of an analyst). The first approach to be discussed is an unsupervised 
computer processing algorithm to detect boundary and interior pixels. The 
interior pixels are classified into unlabeled categories. Then, individual 
boundary pixels and their adjacent pixels are automatically examined to deter- 
mine the categories in which they overlap. The boundary pixels are then 
modeled as a percentage of each adjacent category by using the mean spectral 
values of those categories. The analyst labels samples of the interior 
pixels in order to label the categories. Thus, boundary pixels would be 
automatical ly detennined and allocated to the adjacent categories on the 
basis of the spectral data, and the analyst would be relieved of having to 
allocate boundarie'j to categories. The possible misgivings about this ap- 
proach are as follows: 

1. It has nut been shewn that an automatic ooundary detection algorithr. 
could locate boundaries with sufficient accuracy. In particular, how 
:nanv actual boundaries will be overlooked, how many interior pixels will 
be mistakenly labeled boundary, and what detrimental effects on pro- 
portion estimation will these have? 

'^he first three techniques are suggested in the minutes of the Procedure 1 
Review Conferen.'p held at Texas ASM University on July 13-14, 1978. The min- 
utes ire contained in a letter to R. P. Heydorn from T. C. Minter (dated 
January 18, 1979). Some of the alternatives are from a private communication 
with R. K. terinington resulting from his attendance at the Conference on 
Maximum likelihood Clustering held at Texas AJM University on February 8-9, 

1 979. 
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2. A basic prpniisr of this approach is Ihat an unsupt-rvi sed clustering 
algorithm can partition a Landsat scene into clusters, each of high 
purity (same category). This is still to be demonstrated. Additionally, 
there is no theoretical reason to believe that omission errors would 
balance commission errors so that unbiased proportion estimates would 
resul t. 

3. A boundary pixel may be boundary on more than two categories. Properly 
formulated, though, the modeling could easily be extended to three or 
four categories. 

In the second approach, whenever a boundary pixel is encountered, the analyst 
would consult an unsupervised mul t i Lempora 1 cluster map for determining which 
ficla Lo assign the boundriry pixel. The analyst would then label the bound- 
ary pixel according to the category of the field. This would provide an 
arbiLrary. objective moans for the allocation of boundary pixels to cate- 
gories. The underlying premises are that the clusters from the chosen 
clustering alaorithm would correspond to categories and that boundary pixels 
would be assigned to the most likely category by the algorithm as it assigns 
pixels to categories on the basis of spectral distance (some algorithms may 
also incorporate spatial distance). 

It is env I ..loned I h.it this '.econd approach would be suitable for segment 
inventory i.roredures whifh use unsupervised clustering followed by either 
sLrci ifitrt .irc.il estimation or cluster lalieling with analyst-labeled samples. 
The pro'Pdural flow would be (1) unsijpcrvised clustering of the image data; 
(2) .iidec. hy a cluster map, .inalyst lahelinq of t h? machine-selected pixels; 
and (3) another [irocos'.ing , such as stratified areal estimation, to provide 
the proportion estimates. 

Possible misgivings arc (I) the lack of a , oretical reason to believe that 
the CiU'.tiis will allocate the pixels in an unbiased manner, (2) the possible 
dependence of the allocation on field si^e, (3) the proportions of each crop 
in tile scene, (4) field shape, (5) ttie particular clustering algorithm used, 
(b) the spectral sep.i ratn' I i ty of the caterjories of interest, and (7) the 
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extent to which boundary pixels form their own cluster. It is not clear, 
though, how important these considerations are. This technique has not been 
investigated. It could prove to he an easy, practical technique for allocat- 
ing boundary pixels to categories. 

4 . 2 fW X 1 f4UM-_L Uf.L 1 JC™ 

With these techniques, an analyst, would be asked to label only interior pixels. 
These pixels would b-? utilized as training samples for a maximum likelihood 
classification of all the pixels within the scene. Thus, the boundary pixels 
would have been allocated to the most likely category on the basis of prob- 
abilities computed from ttie labeled pixels. This technique has the same 
advantage as that of clustering; it offers an arbitrary, objective means 
for tiie allocation of boundary pixels to categories. 

from the standpoint of an efficient use of labeled training samples, this 
technique would be besi. suited to proportion estimation of the categories of 
interest by the processing machine. The procedural flow would be for the 
analyst to (1) rereive imagery and spectral aids, (?) label interior pixels 
sampl’d bv some prescribed sch.me, and (3) submit them for classification 
which would provide the estimate. Three variations on how the estimate could 
be produced are as follows; One vaiiation would be to count the number of 
pixels classified info eacli category of interest. A second variation would 
be to Lfmp'.ite t!ie probability for each category of interest for each pixel 
and tliLMi average the probabi 1 i t i es over the scene. This would not provide a 
classification map but may provide more accurate estimation. A third varia- 
tion would be the use of a boundary detection algorithm to separate the 
interior p'i)‘C'l^ f*'hn the boundary ones. Then, pixel counting would be 
used on tne in ucior pixels and probability averaging on the boundaries. 

As can he oasily seen, the number of [<oss ibi 1 i ties here are limited only by 
imaginat ion. 

Ihe possible misgivings are (1) a lad of a theoretical reason for believing 
that tne pixels would be allocated in an unbiased manner, (2) the allocation's 
pos iblr' d^’ptTidence on the proportions of each crop in the scene, (3) the 








■ 

[' 
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dependence on the particular classification algorithm used, and (4) the spec- 
tral separability of the categories of interest. 

4.3 STATISTICAL SAM P LING TECHNIQUE 

With this technique, the analyst would no longer be required to place a single 
category label on a boundary pixel. Instead, he would provide two labels — 
one for each of the categories straddled by the boundary pixel. Each cate- 
gory as labeled by the analyst would then contain both interior pixels and 
boundary pixels. This method of labeling is best suited for use with seg- 
ment inventories obtained by stratified areal estimation. The usual formulas 
have to be modified slightly to accommodate the boundary labels. The deriva- 
tion for the two-category case follows. 

4.3.1 EXAMPLE OF STRATIFIED AREAL ESTIMATION FORMULA MODIFICATION FOR 
BOUNDARY PIXELS 

The usual formula for the proportion estimate from stratified areal estima- 
tion for two categories is the following: 

P = -IW + "_2W 

“ Base n.| Base n^ 

where 

P^j = The estimated proportion for category W, one of the two categories. 

N^. = The number of pixels in the segment classified into class i, 
where i = 1, 2. 

Base = The number of pixels within the entire segment less any unidentifiable 
areas such as clouds. 

n- = The number of labeled samples classified into the it?: class, where 
i = 1, 2. 

= The number of samples labeled W and classified in class i, where 
i = 1, 2. 

The required modification is based on the observation that the n^. pixels in 
the i.'?; class actually consist of some interior pixels n.p and some boundary 
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pixels n.g. Thus, n. = n^p + n^g for i =1,2. The interior pixels n^p are 
further divided into those that are labeled those labeled (j) 

(the second category). These interior pixels n^p^ are wholly within the W 
category. However, the boundary pixels for the ith class are considered only 
partly within the W category; that is, a fraction of the W category. 

Thus, in the formula, n^^ is replaced by n^p^ + 3.|-n^gy. Thus, the modified 
formula becomes: 

p - (oipvj ^ (Ogpyj + ^2'^2BW^ 

w " Base n^ Base n^ 

In this form, the is the average percentage of the boundary pixels classi- 
fied into the \t,i class that should be allocated to the W category. 

4.3.2 DESCRIPTION OF A PRELIMINARY STUDY OF THE TECHNIQUE 

Among the inputs needed to apply this modified formula are the classification 
percentages, the number of labeled test samples classified into each of the 
two categories, and the labels for the test samples as W, (J>, or B. All of 
these are readily available after the processing of a segment. In addition, 
the two 6's are required. The 3's cannot be estimated for each segment on 
a real-time basis, but they can be estimated in advance using blind site 
data. In addition, it is possible to adjust the 3's for use on a per-segment 
basis by making the 3's a function of segment characteristics. Two consider- 
ations for adjusting the 3's for each segment are to (1) select variables 
related to the 3's and (2) quantify the variables for nonblind sites. The 
3's arc thought to be functions of field size, shape, and separability for 
the category of interest. Thus, they may vary from region to region. They 
should, however, be fairly stable from year to year in the same area. This 
allows the use of the blind sites to estimate the 3's for a given area. 

Though it was not attempted in this study, the ratio of the number of 
boundary pixels to the number of pixels labeled is a potential variable that 
is correlated well with the S's. This variable is readily available after 
segment dot labeling. The procedure for determining the 3's for a particular 
area is simply to use ground-truth labels in lieu of analyst-labeled pixels 

A 

and substitute the ground-truth petcentage for the proportion estimate P^. 
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A machine processing (classification) is performed, and the results are input 
to the formula. This provides an equation in two unknowns 0^ and Sg for each 
blind site processed. A constrained regression analysis is performed to 
provide unbiased estimates for 3-| and 02 * assumption in estimating 0-| 
and 02 in this manner is that they are constant for all the blind sites in 
the regression analysis. 

A small study of this technique was conducted. Section 4.3.3 describes 
the data set and the results. The data are presented in the appendix. 
Basically, 26 blind sites in North Dakota and 19 in Kansas were used to 
estimate 3's separately for each state. Because different definitions of 
boundary pixels will result in different 6's, two different definitions were 
considered, one suitable for Procedure 1 and the other suitable for LIST. 

The results indicated that the variances of the proportion estimates of 
small grains in the boundary pixels (0's) are unduly large. In part, this 
is due to the very difficult, i^ not impossible, task of determining the 
"true" proportion of small grains within a set of 209 dots. The difficulty 
is that even with the use of ground truth, there is and must be some sub- 
jectiveness in determining a label for some of the dots. In fact, there 
really is no correct label for the dots which are partially small grains and 
no way to measure directly the percentage of small grains in these partially 
small -grain pixels. Thus, the true proportion cannot be directly measured 
without error. Consequently, the small-grain proportion for the entire 
segment was used. Due to sampling variance, the small-grain ground-truth per- 
centages for the entire segment were sometimes physically irreconcilable 
with the ground-truth percentages for the 209 grid dots. Specifically, 
either the small -grain percentage given by only the interior dots from the 
209 was more than the segment proportion, or the small-grain percentage given 
by the number of interior dots plus all the boundary dots (treated as 100 per- 
cent small grains) in the 209 was less than the segment proportion. Further 
analysis to obtain more reliable estimates has not yet been undertaken. 
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4.3.3 DESCRIPTION OF THE DATA SET AND RESULTS 


A group of analysts was asked to interpret the 209 pixel -grid intersections 
for 45 Phase III blind sites, 26 in Kansas and 19 in North Dakota, to identify 
interior and boundary pixels. The following labels were used: 


Label 

Definition 

A 

Anomalous pixel* 

D 

Nonagricultural area 

P 

Interior pixel 

R 

Misregistered pixel 

X 

Clouds and cloud shadows 

1 

Boundary (small grain and nonsmall grain) 

2 

Boundary (nonsmall grain and nonsmall grain) 


♦Pixel is not representative of most of the other 
pixels within the field; e.g., a mud puddle in a 
wheat field. 


The following blind 
1602 
1604 
1606 
1616 
1619 
1622 
1625 
1635 
1637 
1640 


sites in North 
(Mountrail ) 
(Renvi lie) 
(Ward) 
(Cavalier) 
(Grand Forks) 
(Ramsey) 

(Dunn) 

(Sheridan) 

(Stutsman) 

(Barnes) 


Dakota were used in the study: 
1648 (Bowman) 

1652 (Stark) 

1661 (McIntosh) 

1663 (Richland) 

1899 (Walsh) 

1902 (McKenzie) 

1903 (Mercer) 

1913 (Hettinger) 

1927 (Sargent) 
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The following blind sites in Kansas were used: 


1032 (Wichita) 

1293 

(Meade) 

1033 (Clark) 

1295 

(Osborne) 

1153 (Jewell) 

1297 

(Dickinson) 

1155 (Phillips) 

1885 

(Rice) 

1158 (Washington) 

1340 

(Sumner) 

1166 (Lyon) 

1343 

(Riley) 

1170 (Harper) 

1346 

(Geary) 

1*75 (Sedgwick) 

1851 

(Graham) 

1180 (Cherokee) 

1853 

(Ness) 

1183 (Labette) 

1859 

(Hamilton) 

1279 (Cheyene) 

1861 

(Kearny) 

1285 (Logan) 

1864 

(Stanton) 

1290 (Ford) 

1881 

(Ellsworth) 


Since individual pixels can change with multiple acquisitions from one cate- 
gory to another from acquisition to acquisition because of misregistration 
in the North Dakota blind sites, a reference date was specified; and the 
analysts used three selected additional acquisitions, which included ripe 
and harvested dates when possible, to label the 209 dots. However, in the 
Kansas blind sites, a reference date with the most distinct field boundaries 
was selected in addition to the three other acquisitions used for labeling. 

The corresponding 209-dot ground-truth labels were used to identify the in- 
terior small-grain dots from the nonsmal 1 -grain dots. The classification map 
from Procedure 1 was used to identify the category in which the dots were 
placed. The following quantities were determined: total small-grain dots 

classified as small grains, total small-grain dots classified as nonsmall 
grains, number of dots classified as small grains, number of dots classified 
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as nonsmall grains, total small -grain/nonsmall-grain boundary dots classified 
as small grains, and total small -grain/nonsmall -grain boundary dots classified 
as nonsmall grains. 


To estimate the effect of boundary pixels on the small-grain proportion 
(i.e., to determine the fraction of boundary pixels that should be considered 
small grains), the parameters 3-| and were computed from the following 
equation. 



Base 


'SGI 


+ B 


'B1 


1 n, 


Base 


f"SG2 . - "B2' 

* ®2 ft; 


where 


Base 

^G1 

"SG2 

"1 

"2 

"bi 

"B2 


Ground-truth small-grain proportion. 

All pixels classified as small grains. 

All pixels classified as nonsmall grains. 

The number of pixels within the entire segment less any unidentifiable 
areas such as clouds. 

Total small-grain dots classified as small grains. 

Total small-grain dots classified as nonsmall grains. 

Number of dots classified as small grains. 

Number of dots classified as nonsmall grains. 

Total smal 1-grain/nonsmall-grain dots classified as small grains. 

Total small-grain/nonsmall-grain dots classified as nonsmall grains. 

Fraction of the boundary pixels which are small grains, given that 
the boundary pixels are classified as small grains. 

Fraction of the boundary pixels which are small grains, given that 
the boundary pixels are classified as nonsmall grains. 


To solve for the parameters and B 2 . a linear regression on and X 2 on 
Y was fitted through the origin with the following model: 

Y = + 32^2 
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where 




Two sets of data were analyzed: Procedure 1 data in which A, D, P, and R 
labels were considered interior pixels and the 1 and 2 labels were considered 
boundary pixels; and LIST data, in which D, P, and 2 labels were considered 
interior pixels and the A, R, and 1 labels were considered boundary pixels. 
Pixels labeled X were excluded from the analyses. 


Table 2 presents the regression coefficients along with the corresponding 
standard deviations separately for each labeling method and state. Results 
indicated that the fraction of small grains in the boundary pixels that are 
classified as small grains is 72 percent for Procedure 1 and 49 percent for 
LIST. For Procedure 1 in North Dakota, 72 percent of the boundary pixels 
classified as small grains was small grains, and 19 percent of the boundary 
pixels classified as small grains was small grains. However, in Kansas, 

72 percent of the boundary pixels classified as small grains was small 
grains, and 3 percent of the boundary pixels classified as nonsmall grains 
was small grains. Extremely large standard deviations of the regression 
coefficients '"eflect the variability in the 209-dot sampling and the varia- 
bility due to the number of segments within the state. For example, segment 
1663 in North Dakota (Procedure 1 data) had a ground-truth wheat percentage 
of 51.84. There were 97 interior small-grain dots and five boundaries. The 
random estimate for wheat from the interior pixels was (97/209) x 100 = 46.41 
percent, and the random estimate for wheat from interior and boundary pixels 
was (102/209) x 100 = 48.80 percent. The ground-truth percentage (51.84 per- 
cent) is greater than the percentage of wheat from the random estimate (48.80 
percent) with the interior and the boundary pixels, indicating that the boundary 
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pixels would need to be considered as greater than 100 percent small grains. 
This variability is due to random sampling. 

TABLE 2.- REGRESSION COEFFICIENTS WITH THE CORRESPONDING 
STANDARD DEVIATIONS 


State 

^^1 

Standard 
deviation 
of 3-| 


— 

Standard 
deviation 
of e >2 

R^ 

coefficient 

of 

determination 

Procedure 1 

North Dakota 

0.72656 

0.54149 

0.19814 

0.34859 

0.44986 

Kansas 

0.72404 

0.99781 

0.03050 

0.30608 

0.09065 

LIST 

North Dakota 

0.49496 

0.19963 

0.43250 

0.16310 

0.78673 

Kansas 

0.58680 

0.49548 

0.54625 

0.47139 

0.35887 


4.3.4 PROPOSALS FOR ADDITIONAL ANALYSIS OF THE EXISTING DATA 

One possible way of reducing the variance of the estimated 3's would be to 
incorporate prior knowledge into the estimation process. In this case, 
when the ground-truth proportion is lower than the random sample esti- 
mate for pure small grains then could be used in place of Pgj for 
estimation of the B‘s. Likewise, when P„.p is greater than the random sample 
estimate for pure small grains plus all boundary small-grain P^j, the P^ 
should be used in place of Pg^ in the regression. The resulting regression 
on the truncated proportions will perhaps bias the estimates of the S's. 
However, they will have smaller variances and will dampen the effect of 
sampling anomalies. 

Another interesting analysis that could be performed is to use the Procedure 1 
type 2 dots to obtain a simple random sample proportion for small grains with 
the analyst labels, which include designations of which pixels are boundary 
to perform the regression estimation for the 3‘s. Tl.cse 3's would estimate 
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the analyst's opinion of the amount of small grains in the boundary pixels 
for each category of machine classification. Small-grain proportion esti- 
mates could then be calculated using the estimated g's and would perhaps 
demonstrate a smaller variance on the proportion estimate than the Procedure 1 
estimate although there should be no improvement in the bias. 

Both of these would provide proportion estimates with equivalent bias and 
perhaps lower variances than the Procedure 1 estimates. Superior approaches 
to estimation of the B's exist, but they entail the collection of new data 
from additional follow-on studies. Recommendations are presented in sec- 
tion 5. 

4.3.5 DISCUSSION OF ASSUMPTIONS AMENABLE TO STATISTICAL TESTING 

The statistical sampling technique for handling boundary pixels represents 
an entire class of techniques corresponding to different assumptions about 
6 -| and 62 for two categories. The following table lists five assumptions 
on the S‘s and explains how they correspond to methods sometimes proposed 
for handling boundaries. 


Assumption 


Correspondence 


= ^2 = 

3 ] = 3 ^ “ where 3 is a constant 

^i f 1 ^ 2 ^ where 3^ and are 
constants 


Each boundary pixel is considered as 
containing 50 percent of the category 
of interest regardless of its classifi- 
cation. 

Each boundary pixel is considered as 
containing a fraction fc of the category 
of interest regardless of its classifi- 
cation. 

The boundary pixels classified into the 
\i.h class are considered as containing 
a fraction 3 ^ for the category of in- 
terest. The Traction will be different 
depending on its classification but will 
otherwise remain constant. 
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Assumption 

Correspondence 

"iPW 

"iPW ^ "iP^ 

The fraction B.j of the category of In- 
terest in the boundary pixels which are 
classified Into the category Is 

considered to be the same as the ratio 
of interior pixels for the category of 
interest to all interior pixels classi- 
fied into the Hh class. 

is subjectively determined and 

varies from segment to segment 
and from analysis to analysis 

This is the present method; by this 
method, an analyst is forced to make a 
determination (label) on each boundary 
pixel . 


One of the biggest advantages of the statistical sampling technique is that 
assumptions like these concerning the (?'s can be tested. 

As was discussed earlier, the definition of boundary pixels is critical 
to any method that attempts to handle them. Thus, assumptions which may be 
suitable for some definitions may not be suitable for others. A main point 
to be made, though, is that the statistical sampling method offers a stand- 
ard against which to measure all other methods for handling boundary pixels. 
No matter what definition is chosen for boundary pixels by a candidate 
method, fhe statistical sampling method can use that definition and provide 
a minimum variance unbiased estimate of the average fraction of a boundary 
pixel that contains the category of interest (even conditioned on the classi- 
fication, if one is given). Hence, a best estimate for 6.j is available and 
can be used for hypothesis testing to facilitate the comparison of methods 
for handling boundary pixels. 

4 . 4 MlJL7if.ATECMiY_ imil NG 

As was mentioned earlier in section 2, the analyst tends to be conservative 
in labeling a category of interest, especially for the boundary pixels. In 
the past, analysts were instructed to identify small-grain pixels and label 
any other category as nonsmall grains. The intent of this approach was to 
label the small grains accurately. This produced high accuracies for those 
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pixels identified as small grains. Whenever a pixel was labeled small grains 
in Phase III, the probabilities of being correct were 91.1 percent (ref. 1). 
The analyst did not label all the small grains as small grains. In fact, 
in Phase III, the analyst was able to identify only 78.6 percent of the 
small grains correctly. One of the major sources of this small-grain 
omission error was the preponderance of boundary pixels, which are of course 
partly small grains being identified as nonsmall grains. 

Thus, if the reason for conservative labeling is that the analyst is only 
concerned with one category of interest, then the obvious response is to 
have the analyst identify all categories within the segment. This approach 
was adopted for North Dakota in August of the transition year and later 
extended to all the spring-grain states. Current plans are to use this 
philosophy of identifying all categories in the corn and soybean experiments. 
No results are as yet available on the effectiveness of this change to multi- 
labeling. The cost has increased from the standpoint of interpretation time 
and additional materials required for the analysis. At this time, the 
approach appears to have a significant potential value in reducing the pro- 
portion estimation bias due to boundary pixels and labeling interior pixels. 
Though not related to boundary pixels, an added feature is the study of other 
labeling error sources. Whenever an omission error occurs, it will be known 
as to what category the analyst was confusing with the category of interest. 
This will greatly aid analyst training and feedback. 

5. SUMMARY 

Boundary pixels have been shown to be highly important as sources of pro- 
portion estimation error with Landsat data. The particular definition chosen 
for deciding which pixels are boundary has been found to be critically impor- 
tant to any method of handling the problem. Some of the many methods for 
handling boundary pixels are described in this paper. Two characteristics 
are common to all thv methods. Fir.t. fhere must be a procedure for deciding 
which pixels to call boundary; and serond, there must be an objective rule 
tor processing them. This processing may be counting, performing a regression 
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dnalysif. on their tabulations, or assigning them to classes based on spectral 
measurements. One of the proposed methods, the statistical sampling approach 
to handling boundary pixels, was shown to be a standard by which any of the 
other methods could be compared objectively and provided a general framework 
for viewing the problem. This method also provided a way of producing un- 
biased proportion estimates from labeling techniques that can be applied only 
to interior pixels. 


6. RECOMMENDATIONS 

The following considerations are i"ecommended for use in planning a boundary 

pixel research program: 

a. Conduct an additional analysis of the existing data. 

b. Conduct a study to determine the actual stability of the 8's of the 
statistical sampling method and establish how they might vary with such 
factors as the fraction of boundary pixels. 

c. Test and evaluate several machine algorithms on how well they can handle 
boundary pixels using the statistical sampling method as a standard. 

d. Investigate over larger geographical areas and with a greater number of 
segments the feasibility of having analysts label only Interior pixels 
and simply denote boundary pixels as "boundary." 

e. Evaluate the dependence of the 3's on field size and shape. 
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APPENDIX 


DATA FROM BOUNDARY PIXEL STUDY ON KANSAS AND 
NORTH DAKOTA BLIND SITES 


APPENDIX 


DATA FROM BOUNDARY PIXEL STUDY ON KANSAS AND 
NORTH DAKOTA BLIND SITES 

The tabulation of results for the small boundary pixel study in Kansas and 
North Dakota are presented in tables A-1 throuqh A-6. Tables A-2 and A-3 
contain the raw tabulations of Procedure 1 and LIST for Kansas and A-5 
and A-6 for North Dakota. The transforiDed values amenable for modeling 
using the statistical sampling approach are in table A-1 for Kansas and A-4 
for North Dakota. 


TABLE A-1.- BOUNDARY PIXEL STUDY OF KANSAS BLIND SITES 



Procedure 1 

LIST 

Segment no. 

D 

^2 

Y 

a 

X2 

Y 


1.7308 

1.0318 

4.2006 

4.3270 

4.1272 

7.2295 


0.0000 

4.8520 

1.6170 

0.4902 

1.9408 

-1.7794 


0.4124 

21.5556 

7.6945 

0.2062 

5.8788 

1.8157 

1155 

0.0000 

6.1815 

2.9983 

0.9037 

4.7550 

1.9527 

1158 

0.0000 

0.0000 

0.9189 

0.6484 

0.9598 

1.7358 

1166 

0.4861 

1.9120 

13.9335 

0.4861 

1.9120 

16.3316 

1170 

3.5329 

10.6872 

-0.9777 

9.0846 

2.2265 

4.7522 

1175 

1.1024 

'.3872 

6.0505 

2.7560 

12.4659 

11.3977 

1180 

1.6268 

7.3236 

-0.7518 

3.2536 

6.8928 

9.8254 

1183 

1.1606 

1.8664 

-0.9446 

6.3833 

3.2662 

6.3719 

1279 

0.6506 

1.7136 

-3.5894 

5.8554 

4.2840 

3.9796 

1285 

0.9338 

10.5952 

0.5601 

1.4007 

1.4448 

3.3909 

1290 

0.0000 

1 .0380 

-7.0656 

9.0080 

0.5190 

3.8812 

1293 

0.6973 

0.9126 

-2.4968 

0.6973 

1 .8252 

-2.0405 

1295 

4.6650 

9.1466 

5.6433 

8.2705 

7.2210 

10.5389 

1297 

0.7197 

4.1931 

-26.1042 

0.7197 

3.2613 

-22.5891 

1340 

2.2328 

6.8775 

-1.5304 

5.5320 

7.7945 

13.6204 

1343 

3.6030 

2.3735 

-0.8845 

4.2035 

1 .9028 

-2.5377 

1346 

0.0000 

1.4223 

8.0594 

0.6231 

0.4741 

7.2602 

1851 

2.0884 

3.5992 

4.1042 

1.0442 

1.3497 

6.4981 

1853 

8.7244 

31.1483 

3.2513 

8.2112 

4.6490 

5.9622 

1859 

0.9526 

2.8926 

0.2775 

2.3815 

1 . 9284 

-2.1330 

1861 

0.4919 

4.2471 

0.1832 

1.9676 

2.3595 

4.5303 

1664 

2.2076 

4.1778 

-0.6148 

4.4152 

7.8914 

0.6644 

18P1 

4.5990 

8.0495 

6.5900 

2.5550 

5.6320 

C.2845 

1885 

5.5627 

6.4120 

14.4641 

8.0912 

6.8700 

13.6435 


, >= 0.60633® 

‘r2 


•■x 

*1 

„ * 0.22873*^ 
*2 


The correlation coefficient is significantly different from zero at 
the 1-percent level. 

The correlation coefficient is rot significantly different from zero. 
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TABLE A-2.- PROCEDURE 1 STUDY OF KANSAS BLIND SITES' 


Segment 

P 

no. 

BC 



^2 

Base 

"l 

"2 

"SGl 

"SQ2 

14038 

22673 

88 

120(1) 

69 

9(1) 

21923 

22702 

7 

199(3) 

1 

14 

22558 

22793 

5 

202 

3 

31(2) 

18153 

19379 

7 

197(3) 

6 

9 

19603 

22693 

21 

180(2) 

17 

17 

21585 

22688 

10 

199 

5 

12 

8492 

22700 

124 

84(1) 

no 

19 

17783 

22660 

39 

170 

36 

39 

17924 

22732 

26 

183 

17 

30 

18137 

22596 

34 

172(3) 

18 

12 

15479 

22724 

49 

159(1) 

45 

11 

1 7978 

22760 

45 

164 

19 

18 

11681 

21433 

101 

105(3) 

94 

12 

10182 

12127 

23 

184(2) 

13 

13(2) 

15225 

22752 

68 

139(2) 

53 

23(2) 

20392 

22677 

14 

193(1) 

12 

42(1) 

16229 

22688 

51 

156(2) 

47 

70 

16712 

21160 

35 

166(8) 

12 

5 

21858 

22707 

6 

203(1) 

4 

2 

18076 

22698 

39 

177(2) 

18 

20 

14780 

22704 

68 

140(1 ) 

36 

15 

17209 

22^32 

51 

157(1) 

32 

29 

12427 

22701 

92 

116(1) 

57 

15 

12520 

18860 

61 

143(5) 

50 

19 

18015 

22644 

40 

168(1 ) 

17 

18 

11334 

22701 

99 

109(1 ) 

67 

13 


’Figures in parentheses represent the number of thrnsholdod pixels. 


"bi 

'’82 

4 

2 

0 

10 

2 

44 

0 

13 



B, * 0.72404 
8, • 0.03050 


0.99781 

0.30608 
























TABLE A-3.- LIST BOUNDARY PIXEL STUDY OF KANSAS BLIND SITES' 


Segtnent 

no. 



A 

p 

^BC 


N2 

Base 

"l 

38.7 

8635 

14038 

22673 

88 

8.9 

779 

21923 

22702 

7 

23.5 

235 

22558 

22793 

5 

12.7 

1226 

18153 

19379 

7 

20.1 

3090 

19603 

22693 

21 

22.1 

1103 

21585 

22688 

10 

63.0 

14208 

8492 

22700 

124 

43.9 

4872 

1 7788 

22660 

39 

26.0 

4808 

1 7924 

22732 

26 

15.1 

4459 

18137 

22596 

34 

30.4 

7245 

15479 

227?*, 

49 

18. 1 

4782 

17978 

22760 

45 

41.5 

9752 

11681 

21433 

101 

12.5 

1945 

10182 

12127 

23 

42.5 

7527 

15225 

22752 

68 

2.1 

2285 

20392 

22677 

14 

56.8 

6459 

16229 

22688 

51 

8.7 

4448 

16712 

21160 

35 

11.5 

849 

21858 

22707 

6 

22.5 

4622 

18076 

22698 

39 

28.7 

7924 

14780 

22704 

68 

29.5 

5523 

1 7209 

22732 

51 

35.3 

10274 

12427 

22701 

92 

35.8 

6350 

1 2520 

18860 

61 

23.8 

4629 

18015 

22644 

40 

54.3 

11367 

11 334 

22701 

99 



^Figures in parentheses represent the number of thresholded pixels. 

0, •= 0.58680 Sfl * 0.49548 

1 


02 • 0.54625 


0.47139 


























TABLE A-3.- LIST BOUNDARY PIXEL STUDY OF KANSAS BLIND SITES*’ 


Segment 

no. 

^BC 

«1 

^-2 

Base 

"1 

"2 

"SGI 

"SG2 

"bi 

"B2 

1032 

38.7 

8635 

14038 

22673 

88 

120(1) 

62 

9 

10 

8 

1033 

8.9 

779 

21923 

22702 

7 

199(3) 

1 

21 

1 

4 

1153 

23.5 

235 

22558 

22793 

5 

202 

3 

43 

1 

12 

1155 

12.7 

1226 

18153 

19379 

7 

197(5) 

4 

15 

1 

10 

1158 

20.1 

3090 

1960.3 

22693 

21 

180(2) 

15 

18 

1 


1166 

22.1 

1103 

21585 

22688 

10 

199 

4 

8 

1 


1170 

63.0 


8492 

22700 

124 

84(1) 

96 

22 

18 


1175 

43.9 


17788 

22660 

39 

170 

33 

31 

5 

27 

1180 

26.0 


17924 

22732 

26 

183 

13 

13 

4 

16 

1133 

15.1 


18137 

22596 

34 

172(3) 

7 

10 

11 

mm 

1279 

30.4 

7245 

15479 

22724 

49 

159(1) 

36 

7 

9 


1285 

18.1 

4782 

1 7978 

22760 

45 

164 

15 

16 

3 

3 

1290 

41.5 

9752 

11681 

2143.3 

101 

105(3) 

72 

10 

20 

1 

1293 

12.5 

1945 

10182 

12127 

23 

184(2) 

13 

12 

1 

4 

1295 

42.5 

7527 

15225 

22752 

68 

139(2) 

37 

29 

17 

15 

1297 

2.1 

2285 

20392 

22677 

14 

193(1) 

11 

36 

1 

7 

1340 

56.8 

6459 

16229 

22688 

51 

156(2) 

33 

54 

10 

17 

1343 

8.7 

4448 

16712 

21160 

35 

166(8) 

10 

11 

7 

4 

1346 

11.5 

849 

21858 

22707 

6 

203(1 ) 

3 

5 


1 

1851 

22.5 

4622 

18076 

22698 

39 

177(2) 

16 

17 


3 

1853 

28.7 

7924 

14780 

22704 

68 

140(1) 

28 

18 


10 

1859 

29.5 

5523 

17209 

22732 

51 

157(1' 

32 

34 


4 

1861 

35.3 

10274 

12427 

22701 

92 

116(1) 

52 

11 


5 

1864 

35.8 

6350 

12520 

18860 

61 

143(5) 

46 

21 

3 

17 

1881 

23.8 

4629 

18015 

22644 

40 

168(1) 

21 

27 

5 

12 

1885 

54.3 

11367 

11334 

22701 

99 

109(1 ) 

65 

17 

16 

15 


^Figures 1n parpnthp^es represent the number of thresholded pixels. 


8, ' 0.58680 s. • 0.49548 

I 8, 



8j - 0.54625 


0.47139 






















TABLE A-4.- BOUNDARY PIXEL STUDY OF NORTH DAKOTA BUND SITES 



Procedure 1 

LIST 

Segment no. 

1 

^1 


Y 

^1 


Y 

1602 

11.6679 

11.7234 

13.2698 

19.7847 

9.9198 

10.9025 

1604 

3.7890 

12.0624 

3.2374 

11 .3670 

27.1404 

18.2598 

1606 

5.0643 

7.6320 

-0.4161 

5.0643 

8.1090 

0.1466 

1616 

3.1044 

5.3364 

3.0688 

P.7958 

11.1175 

10.3936 

1619 

1.8248 

1.0210 

0.7811 

3.1934 

3.0630 

3.6269 

1622 

1.8040 

1 . 0308 

2.6032 

8.5690 

4.1232 

7.4998 

1625 

0.0000 

2.8566 

-1.7298 

0.0000 

4.7610 

-1.5924 

1635 

0.0000 

0.5266 

-3.4571 

0.0000 

0.5266 

-4.6316 

1637 

0.0000 

0.4727 

-3.5637 

0.0000 

0.4727 

-2.6183 

1640 

1 .1018 

1.3014 

-2.3832 

4.9581 

6.9408 

4.9781 

1648 

2.5242 

5.0670 

5.6236 

3.7863 

8.1072 

6.6370 

1652 

3.2697 

6.3817 

-0.0360 

7.4736 

15.2179 

7.1847 

1661 

0.9490 

4.4838 

2.0305 

0.4745 

0.0000 

-2.4059 

1663 

0.9918 

1 .4196 

4.5095 

0.4959 

1.4196 

3.5858 

1889 

1.3968 

0.0000 

0.7601 

12.1056 

3.0402 

8.3741 

1902 

0.0000 

3.8112 

2.9232 

0.0000 

0.4764 

1 .0176 

1903 

1 .6761 

5.6280 

5.9761 

8.9392 

8.9110 

12.5011 

1913 

0.5604 

2.7060 

11.8000 

2.8020 

5.9532 


1927 

0.0000 

0.9390 

-3.4258 

0.5146 

2.8170 


r 

V „ = 0.77586^ 

'r2 


^1 

y = 0.58203^ 


^The cnrrela Lion coefficient is significantly different from zero at the 
1 -percent level. 
















TABLE A-5.- PROCEDURE 1 BOUNDARY PIXEL STUDY OF NORTH DAKOTA^ 


Segment 

no. 

^GT 

Nl 

N2 

Base 

"l 

"2 

"SGI 

'’SG2 

"81 

"B2 

1602 

38.24 

6534 

6092 

12626 

102 

107 

43 

7 

23 

26 

1604 

52.42 

5930 

16784 

22714 

62 

147 

44 

61 

9 

24 

1606 

32.94 

3937 

10637 

14574 

48 

153(8) 

33 

31 

9 

16 

1616 

66.76 

11403 

11320 

22728 

97 

112 

87 

42 

6 

12 

1619 

52.69 

11848 

10931 

22779 

114 

94(1) 

97 

15 

4 

2 

1622 

5G.28 

12271 

10400 

22671 

120 

89 

92 

12 

4 

2 

1625 

21.54 

1725 

20890 

22615 

14 

194(1) 

13 

34 

0 

6 

1635 

16.03 

2697 

6808 

9505 

70 

136(3) 

13 

27 

0 

1 

1637 

35.77 

3992 

14552 

18544 

36 

166(7) 

31 

44 

0 

1 

1640 

52. 1C 

10872 

11810 

22682 

87 

120(2) 

80 

24 

2 

3 

1648 

20.29 

5453 

17282 

22735 

57 

150(2) 

18 

14 

6 

10 

1652 

30.62 

5090 

17609 

22699 

48 

158(3) 

32 

32 

7 

13 

1661 

40.82 

6439 

13330 

19819 

69 

135(5) 

45 

35 

2 

9 

1663 

51.84 

7759 

14915 

22674 

69 

139(1) 

63 

34 

2 

3 

■n 

59.33 

15205 

7469 

22674 

144 

65 

116 

9 

3 

0 

mm 

8.64 

315 

22445 

22760 

1 

207(1) 

0 

12 

0 

8 

mm 

17.35 

2792 

19923 

22715 

22 

187 

17 

4 

3 

12 

1913 

29.89 

2675 

20051 

22726 

42 

163(4) 

24 

21 

2 

5 

1927 

31.36 

7207 

15378 

22585 

62 

154(2) 

53 

16 

0 

2 


^Figures in parentheses represent the number of thresholded pixels. 
0, • 0.72656 So ' 0.54149 

B, • 0.19814 So • 0.34859 

2 02 


y- 


f 





















TABLE A-6.- LIST BOUNDARY PIXEL STUDY OF NORTH DAKOTA'’ 


Segment 

no. 

^GT 


38.24 

■■ 

52.42 

1606 

32.94 

1616 

66.76 

1619 

52.69 

1622 

50.28 

1625 

21.54 

1635 

16.03 

1637 

35.77 

1640 

52 10 

1648 

20.29 

1652 

30.62 

1661 

40.82 

1663 

51.84 

mam 

59.33 


8.64 

■■ 

17.35 

1913 

29.89 

1927 

31.36 



3937 

11408 

11848 

12271 

1725 

2697 

3992 

10872 

5453 

5090 

6489 

7759 

15205 

315 

2792 

2675 

7207 


11320 
10931 
10400 
20890 
6808 
14552 
11810 
17282 
17609 
1 3330 
14915 
7469 
22445 
19923 
20051 
15j78 


12626 

22714 
14574 
22728 
22779 
22671 
22615 

9505 

18544 

22682 

22735 

22699 

19819 

22674 

22674 

22760 

22715 
22726 
22585 


"2 

"SGl 

107 

45 

147 

31 

153(8) 

32 

112 

78 

94(1) 

93 

89 

80 

194(1) 

11 

136(3) 

12 

166(7) 

31 

120(2) 

69 

150(2) 

13 

158(3) 

26 

135(5) 

47 

139(1) 

62 

65 

104 

207(1) 

0 

187 

7 

163(4) 

19 

145(2) 

52 



’Figures in parentheses represent the number of thresholded pixels. 


e, = 0.49496 


= 0.43250 


s,, = 0.19963 
s., = 0.16310 
















